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Abstract —In this paper, interference alignment (IA) is de¬ 
signed for secure multiuser multiple-input multiple-output sys¬ 
tems in the presence of an eavesdropper. The proposed IA 
technique designs the transmit precoding and receiving subspace 
matrices to minimize both the total inter-main-link interference 
and the wiretapped signals. With perfect channel state infor¬ 
mation is known at the transmitters (Txs), the cost function 
of the optimization problem is alternatively minimized over 
the precoding matrices at the Txs and the receiving subspace 
matrices at the receivers (Rxs) and the eavesdropper. The feasible 
condition, the proof of convergence of the proposed IA approach 
are provided. Numerical results reveal a significant improvement 
in terms of average secrecy sum rate of our IA algorithm as 
compared to the conventional IA design. 

I. Introduction 

Due to the broadcasting nature of the wireless multiuser 
multiple-input multiple-output (MIMO) communication net¬ 
works, the data transmission among legitimate users is sus¬ 
ceptible to be wiretapped by nearby eavesdroppers. Various 
signal processing techniques in the physical layer have been 
developed to enhance secrecy against the eavesdropper [1]- 
[5]. However, all the solutions in [1]—[5] mainly focuses on the 
optimization problems of secrecy sum rate (SSR) in the broad¬ 
cast channels in which only one base station communicates 
with multiple users in the network. The secure communication 
for multiple transmitter-receiver (Tx-Rx) pair transmissions 
has not been studied. It should be emphasized that the sum 
rate maximization of the MIMO interference channels with 
multiple user pairs is mathematically challenging even for the 
scenarios without taking secure transmission into considera¬ 
tion. Interference alignment (IA) has recently been exploited 
as a potential technique to increase the secrecy in multiuser 
MIMO systems [6], [7]. The authors in [6] has proven that it 
is possible for each user in the network to achieve a nonzero 
secure degrees of freedom (DoF) when using an IA scheme 
along with precoding matrix at each Tx. It has been shown 
in [7] that the secure transmission is also feasible where the 
number of antennas at legitimate Tx ans Rx is greater than 
those of the eavesdropper. However, a detailed IA design for 
secure multisuer MIMO communication networks has not been 
explored in the literature. 

In this paper, we study on the secure communication for 
multiuser MIMO systems in presence of an eavesdropper. 
Inspired by the standard IA technique in [8]—[ 10], we de¬ 


velop the IA scheme, namely the wiretapped-signal-leakage- 
minimization (WSLM) IA technique, for the secrecy trans¬ 
mission in multiuser MIMO systems. By forcing the receiving 
signals at the eavesdropper into reduced-dimensional receiving 
subspace, and then minimizing the total signal power, the 
desired signal will become harder to be recoverd at the 
eavesdropper. The standard IA technique has been well studied 
in MIMO systems without the eavesdropper [8]—[10]. The 
authors in [8] have proposed a scheme which make the Rxs 
interference-free by forcing interfering signals into a reduced¬ 
dimensional subspace. However, closed-form solutions of the 
precoders with more than 3 users appear to be critical. There¬ 
fore, iterative techniques have been investigated [9], [10], in 
which the precoding matrices at Txs and the interference sub¬ 
space matrices at Rxs are updated in every iteration. Compared 
with the scheme in [8], the iterative algorithm can efficiently 
obtain the optimal precoding matrices and interference sub¬ 
spaces for the system with more than 3 Tx-Rx pairs. Motivated 
by these works, we adopt the iterative approach to design IA 
for the secure multiuser MIMO systems. We also derived the 
proof of convergence, the feasible condition and the analysis 
on system parameters for the best IA design. Simulation results 
show that the proposed IA design outperforms the conventional 
approach [10] in terms of average SSR. 

The rest of this paper is organized as follows. Section II 
introduces the multiuser MIMO system in the presence of an 
eavesdropper. In Section III, we propose the IA algorithm. The 
numerical results are provided in Section IV. Finally, Section 
V presents concluding remarks. 

Notation: X and x are respectively denoted as a complex 
matrix and a vector. X 1 and X H are the transposition and 
conjugate transposition of the complex matrix X, respectively. 
/ and 0 are respectively identity and zero matrices with the 
appropriate dimensions, tr(.), rank(.) and E(.) are the trace, 
rank and expectation operators, respectively. ||JT||j? is the 
Frobenius norm, x ~ CAf(x. R x ) means that x is a complex 
Gaussian random vector with means x and covariance R x . 
TUX) and M(X ) are respectively the range and nullity of 
matrix X. X and dim],T} is denoted as the orthogonal 
subspace and the dimension of subspace X , respectively. 

II. System Model 

Consider a MIMO interference channel with K Tx-Rx 
pairs and an eavesdropper as depicted in Fig. 1. Without 
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Fig. 1. A system model of secure multiuser MIMO communication network. 


loss of generality, we assume that the eavesdropper is the 
(K + l)-th Rx with the N re antennas, all Txs and Rxs are 
equipped with the same number of antennas N t and N r , 
respectively, and communicate with the same number of data 
stream d. The received signal y k £ C JV, ’ xl at the k -th Rx for 
k £ K. = {1, K) can be given as 

K 

y k = YH ktl F t s t +n k (1) 

where H k ,i £ C NrXNt is the static flat-fading MIMO channels 
matrix from the i- th Tx to the fc-th Rx, F k £ C NtXd is the 
precoding matrix applied on the signal vector s k £ C dxl and 
rife ~ CN (0, cr^/jv r ) is a complex Gaussian noise vector. 
Assuming that EjsfcSjf} = Id and the global channel state 
information (CSI) is available, the channel capacity at the fc-th 
legitimate Rx, for k £ K,, can be calculated as [11] 

TZ k = log 2 \l Nr + H ktk F k Fj?Hj d k R~^ | (2) 


where [a] + = max(a,0). 

III. WSLM for Secrecy Multiuser MIMO 
Communication via IA 

The aim of WSLM IA design is to find the precoding 
matrices {F k } k=1 which aligns all interference signals at the 
legitimate Rxs into reduced subspaces and all signals at the 
eavesdropper into a subspace. Mathematically, the precoder 
matrices {F k} k ^ and the receiving signal subspace matrix 
{Wk} k +i (at the legitimate Rxs and the eavesdropper) satisfy 
the three following conditions 

rank(Wf H k)k F k ) = d (6) 

W»H k , e F e = 0;W^M€/C (7) 

Wk + 1 HK+i,eFi = 0 ; W € /C. (8) 

Conditions (6) guarantees that the intended signal at the fc-th 

Rx achieves d degree of freedom when H k>k is full rank. The 
condition (7) ensures no inter-main link interference (IMLI), 
i.e., no interference among the legitimate users, at the fc-th 
Rx. Condition (8) makes certain that there is no leakage signal 
wiretapped at the eavesdropper or the (K + l)-th Rx. 

To remove the IMLIs at the intended fc-th Rx, we align 
the interference signal H k ^Fe into the interference receiving 
subspace W,. 1 which is spanned by the orthonormal basis 
matrix U k . Hence, the total interference leakage inside the 
receiving subspace is defined by 

M{F k }« =1 ,{U k }% =1 ) (9) 

K K 

= Y Y \\H k ,eF e -U k U”H k/ F e \\l 

fc=l 


To reduce the wiretapped information, we align all the 
receiving signals at the eavesdropper into a receiving subspace 
Wk -)_i spanned by the orthonormal basis matrix Wk- i-i where 
their total signal power is minimized. Mathematically, we aim 
to minimize 


K 

where R Zk = ^ H k ^F t F^H ki + a\I jv r is the cor- 

i=i,e^ k 

relation matrix of interference and noise in Eq. (1). The 
eavesdropper or the (K + l)-th Rx wiretaps the data signal 
from the k -th Tx-Rx pair and the information leakage rate of 
this wiretap channel can be given by 


ni e) = io g2 


[ N r , 


H 


F T?H ztH n —1 

k* k ^K+l,k^e,k 


(3) 


32{{Fk} k =\,U k+ t) (10) 

K 

= Y \\H K +i,tF t ~ U k+1 U” +1 H k+u Fz\\ 2 f 

fe=l 

where Uk+i is the orthonormal basis of the subspace W^ +1 . 

Finally, the IA problem can be expressed as the following 
joint optimization problem 


K 

where R e , k = ]T H K+ xjF e Ff H% +11 + o\I Nr<i is the 

interference plus noise correlation matrix at the eavesdropper. 
Now, the secrecy rate for the fc-th Tx-Rx pair can be defined 
as 

F s , k = [n k - n[ e) } + - (4) 


min 

{F k }£ = iAu k If+i 1 


s.t. 


where 


3({F k } k =i> {U k} k =i) (11a) 

F*F k = ^I- Wk£lC (lib) 

U k U k = I;Wk £ /C U {K + 1} (lie) 


Finally, the SSR of multiuser MIMO system can be expressed 
as 

K K 

n s = Y 7Zs ’ k = YV lk - n k ) ] + ’ (5) 

k= 1 k= 1 


3({F k }% =1 ,{U k }™) (12) 

= 3i ({F k }X =1 ,{U k }« =1 )+M{F k }£=i,U K+1 ) 

We now derive the solution of problem (11) via alternating 
minimization [10] in the following. 






















Transmit precoder selection: When {U k } k +^ are fixed, 
for finding Ft, (11) reduces to the optimization problem 


min 

Ft 


S.t. 


tr 


K+l 

F? [ e hU 1 


Ff Ft = ~jli W £ K 


U k U*)H k ^j Ft 

(13a) 

(13b) 


The solution to the above optimization problem is given by 

[ 12 ] 


K +i j 

E Hlt{I-U k UX)H k A (14) 

k=l,k^l I 

where Cmin{-^} ' s a matrix whose columns are the d eigen¬ 
vectors corresponding to the d smallest eigenvalues of X. 

Receiver interference subspace selection: Secondly, when 
Ft is fixed, from (11), the orthonormal basis U k can be 
selected by solving the following optimization problem 



Algorithm 1 : Proposed WSLM IA Algorithm for Secure 
Multiuser MIMO Systems 


satisfied 


then select {U y k } k Jf from (16) 


l: Inputs: d,H k ,t, Vfc £ /C U {K + 1}, W £ /C, k = 0, k, 

where k is the iteration index; 

2: Initial variables: random matrix 

jp(0)H p(0) _ /Ft r 

r k c k ~ V d 1 
and (18); 

3: Evaluate the objective function 3({F k 
from (11); 

4: while K < Kmax do 

5: For fixed {U k } k ^, select {F\f +1 ^} k Li from (14); 


ifc= l'\ u k Sk=l ) 


ffc=l 

For fixed {F k +1 }J? =1 , select {(/), 


} k Jf from (16) 

and (18); 

Evaluate the objective function 

3({F ( k +1) }Li,{U[ K+1) Ifjj 1 ) from (11); 
k = k + 1; 

Repeat steps 5-8 until convergence. 

end while 


A. Convergence analysis for the WSLM IA design 


max 

u k 

1 

ul 

s.t. 

ul! 

u k 


K 




Hence, U k can be given by [12] 

K 


u k = C, 


N r -d 

max 


E H k , e FtF”H» 


(15b) 


(16) 


In the K-th iteration, after step 5 and 6 in Algorithm 1, it 
can be verified that 

3({iE +1) }f =1 , {U[ K+1) Ifjj 1 ) < 3 ({^ k) }Li, j 1 ) 

(19) 

That is, the cost function j) in (11) is reduced monotonically 
over iteration. In addtion, the cost function is bounded below 
by zero. Thus, the convergence of the first proposed IA 
algorithm is guaranteed. 

B. Feasible condition for the proposed WSLM IA design 


where £max d {X } is a matrix whose columns are the (N r — d) 
dominant eigenvectors of matrix X. 

To deal with wiretapped signals, we assume that Wk +i 
is a N re x n matrix where n < N Te . Since Wk+ i is an 
orthonomal basis matrix, rank (f/V k+i) = rank(W^ + i) = 
n. To satisfy condition (8), Hx+i,eFt G C Nr ‘ xd must lie 
in the null space of Wf i+1 . Hence, the existence condition 
of Wk +1 is N Te — n > d. Thus, we consider the case that 
n = N Te — d for the dimension of the null space of W^ +1 
to be smallest, to make the eavesdropper harder to recover 
signal. Then, the matrix size of Uk+i must be N Te x d. From 
(11), the orthonormal basis Uk+i can be found by solving the 
following optimization problem 


max tr{[/f +1 [E H K+u F e F?.H% +i e J U K +i 

Uk + 1 ' ’ 


1 =1 


s.t. Uk+iUk+i=I( 1 - 
The solution to (17) can be found as follows [12] 


K 


u K+ i = c; 


d 

max 


E^+t^f^+t,/, 


t=i 


(17a) 

(17b) 


(18) 


The step-by-step iterative algorithm for the WSFM method 
can be shown in Algorithm 1. 


The feasibility of the IA scheme is well studied in [13] but 
not in the security context and, thus, in this subsection, we 
derive the feasible condition for the proposed IA approach. It 
has been shown in [13] that the proper systems surely render 
the feasibility of the IA problems. The MIMO interference 
channels are known as proper systems if the number of 
variables N v is greater than or equal to the number of equations 
N eq in the IA conditions. The total number of equations N eq 
is directly given from (7) and (8) as follows 

N eq = K(K - 1 )d 2 + K(N re - d)d. (20) 

The number of variables to be designed for the precoder Ft, 
subspace matrices W k and Wk+i equals d(N t —d), d(N r — d) 
and d(N re — d), respectively [13]. Finally, the total number of 
variables in the proposed IA conditions (7) and (8) can be 
given as 

N v = Kd (N t + N r - 2d) + d(N re - d). (21) 

Denote (N t x N r , N re , d) K as the system where there 
are K Tx-Rx pairs and one eavesdropper; the Tx, Rx and 
eavesdropper are respectively equipped with N t , N r , N Te 
antennas; and d data streams are sent between each Tx-Rx pair. 
By comparing N eq in (20) and N„ in (21), the IA scheme in 
Algorithm 1 for the (N t x N r , N re , d) K system with secure 
communication is feasible if N v > N eq , i.e., 

K(N t + N r ) - (K 2 + 1 )d> N re (.K - 1). (22) 




Remark 1. Condition (22) implies that for satisfying the 
WSLM IA feasible conditions, the number of antennas at 

, , . . , , K(N t +N r )-(K 2 +l)d 

the eavesdropper is restricted by N re < - 1 -— 


K-l 


Since the number of antennas at the eavesdropper cannot be 
controlled in the proposed system, the secrecy can be enhanced 
when increasing the number of transceiver antennas N t and 

N r . 


C. Relationship between information-leakage-rate and the 
number of antennas at the eavesdropper 

In the WSLM IA design, we aim to minimize the total 
wiretapped signal power in a reduced dimension subspace 
and, thus, the sum information-leakage-rate (SILR), defined as 

K 

TZ ^, can be decreased. However, increasing N re means 

fc=i 

increasing the dimension of the decoding space at the eaves¬ 
dropper and, therefore, results in an increase in the SILR. 
Hence, one can ask which is the relationship between N re 
and the SILR in the WSLM IA design. We investigate the 
answer of this question in this subsection. From (7), we have 

7 Z(F t ) C K{Uj[H kt t),V£?k,£elC (23) 

K(F k ) C Af(U^H kti ),W£^kJelC. (24) 


From (8), the eavesdropper can overhear the signal of the fc-th 
Tx if the subspace spanned by F k lies in the subspace spanned 
by W% +1 H k+ i, fc , i.e, 

K(F k ) C n(W% +1 H K+lik ),VkeK. (25) 

It is noted that 


dim {n(U%H k , t )} 


dim{N (U%H k/ )} 
dim {7 l(W$ +1 H K+ltk )} 


rank ( U k H k ,e ) 

rank (U k ) = N r — d, (26) 

M - dim{7?. 

N t -N r + d, (27) 

rank {W^ +1 H K +i,k) 
rank {W K + i) 

N re - d. (28) 


From (24) and (25), it is implied that the eavesdropper can 
recover the signal of the fc-th Tx easier if the dimension of 
space for recovering this desired signal at the eavesdropper is 
larger than the dimension of the space where this desired signal 
belongs to. Mathematically, the SILR will start increasing 
when 

dim {H(W% +1 H K +i,k)} > dim{W(C7fiT M )},(29) 
44 N re > N t ~N r + 2d. (30) 


Remark 2. It can be seen from (26) and (27) that 
if N r increases, dim{7£ (U k H k } increases and 
dim{A/” (U k H k ,e) } decreases, which makes Fi and U k 
have the larger dimension of the search space and easier 
to reach the optimal solution, for all £,k £ KL. Hence, the 
main-link rate is increased when N r increases. In addition, it 
can be seen from (30) that when N r increases, (N t —N r + 2d) 
decreases, which means the number of antennas required for 
keeping a low SILR at the eavesdropper becomes smaller. 
Therefore, there should be a trade-off between maximizing 
main-link-rate and the SILR for different values of N r . 


IV. Simulation Results 

In this section, we evaluate the SSR performance of our 
proposed IA design through some numerical results in compar¬ 
ison with the conventional IA algorithm [10]. In simulations, 
noise variances are normalized a k = a 2 = 1. The Rayleigh 
fading channel coefficients are generated from the complex 
Gaussian distribution CAf(0, 1). We define signal-to-noise-ratio 
SNR = All the numerical results are averaged over the 200 
channel realizations. 

First, we investigate the convergence characteristic of the 
WSLM IA algorithm. We run the simulation for a (9 x 9,6, 3) 3 
system with a random channel realization. Note that this 
system satisfies the feasible conditions (22). The evolution of 
the cost functions (SSR) over iterations is illustrated in Fig. 
2. As can been seen from Fig. 2, the cost function reduces 
monotonically over iteration and approaches to zero. The 
proposed IA algorithm takes around 25 iterations to converge 
to 10 -9 . 

Now, we compare the average SSR of our proposed IA 
algorithm with that of the conventional IA method [10] which 
does not involve the secrecy constrains. We consider two sce¬ 
narios of (9 x 9, 6,3) 3 and (9 x 9, 9,3) 3 systems. It should be 
noted that the WSLM IA feasible condition is satisfied in these 
two systems. It has been revealed from Fig. 3 that the WSLM 
IA scheme significantly improve the SSR in comparison with 
the conventional IA design when N re is high. That is because 
when N re is small, the SILR is small and hence, the difference 
between our IA design and the conventional IA one is not 
noticeable. 

Finally, we investigate impacts of the number of receive 
antennas N re on the sum main-link-rate (SMLR) and the 
SILR. Fig. 4 illustrates the average SMLR, SILM while 
Fig. 5 plots the SSR of the proposed IA design for the 
(20 x 9,7V re ,3) 3 and (20 x 20,iV T . e ,3) 3 systems with N re = 
3,..., 25 at SNR = 30 dB. As such parameters, the WSLM 
IA feasible condition is satisfied in these two systems. It 
should be noted from (30) that the SILR will start increasing 
if N re > 6 in the (20 x 20, N re , 3) 3 system and N re > 17 
in the (20 x 9,N re ,3 ) 3 system. These analytical results have 
been confirmed by the numerical simulation in Fig. 4. On the 
other hand. Figs. 4 and 5 show the trade-off between the 



Fig. 2. The convergence behavior of the proposed IA algorithm. 






V. Conclusion 



In this work, we design an IA technique for secure mul¬ 
tiuser MIMO system with the presence of an eavesdropper. 
When the CSI of main links and wiretapped links are known 
at Txs, the transmit precoders and the receiving subspace ma¬ 
trices are selected respectively in every iteration to minimize 
both IMLIs and wiretapped signals. Our analysis provides 
key insights into the proposed IA design over the proof of 
convergence, feasible condition and relationship between SILR 
and the number of antennas at Txs anf the eavesdropper. 
The numerical results show that the proposed IA algorithm 
obtained a significant SSR improvement in compared with the 
conventional IA method. 


Fig. 3. The average SSR versus SNR for (9 X 9, 6, 3) 3 and (9 X 9, 9, 3) 3 
systems. 



N (Number of antennas at the eavesdropper) 


Fig. 4. The average SILR and SMLR versus N re for (20 x 9, N re , 3) 3 and 
(20 X 20, N re ,3) 3 systems. 



Fig. 5. The average SSR versus N re for (20 x 9, 7V re , 3) 3 and 
(20 X 20, Nr e , 3) 3 systems. 
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SMLR and SILR for different number of receive antennas 
{N r ,N re }. When N r is large (e.g., N r = 20) the SMLR 
is higher. However, the number of antennas N re required for 
keeping a low SILR at the eavesdropper is smaller (see (30)). 
Thus, the SSR for the case N r = 20 is high at first when 
N r<t is small but then it quickly reduces as N re increases in 
comparison with the case of N r = 9. 














